#include "functions.h"
using namespace std;

void display(double** A, int n)
{
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
        {
            cout << A[i][j] << " ";
        }
        cout << endl;
    }
}

int main()
{
    int n = 3;
    double** A = new double*[n];
    for(int i = 0; i < n; i++)
        A[i] = new double[n]();
    A[0][0] = -4.0;
    A[0][1] = -3.0;
    A[0][2] = -7.0;
    A[1][0] = A[1][2] = A[2][1] = 2.0;
    A[1][1] = 3.0;
    A[2][0] = 4.0;
    A[2][2] = 7.0;
    double** Q = new double*[n];
    for(int i = 0; i < n; i++)
        Q[i] = new double[n]();
    Hessenberg_Householder_Resolution(A, Q, n);
    display(A,n);
    display(Q,n);
    system("pause");
    return 0;
};